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Abstract 


A ring partially filled with a viscous fluid has been analyzed 
as a nutation damper for a spinning satellite. The fluid has been 
modelled as a rigid sl\ag of finite length moving in a tube and resisted 
by a linear viscous force. It is shown that there are two distinct 
modes of motion, called the spin synchronous mode and the nutation 
synchronous mode . Time constants for each mode are obtained for 
both the symmetric and asymmetric satellite. The effects of a stop 
in the tube and an offset of the ring from the spin axis are also 
investigated. An analysis of test restiLts is also given including 
a determination of the effect of gravity on the time constants in 


the two modes . 



Contents 


Abstract 

Illustrations 

Synibols 

1. Introduction 

2. Statement of the Prohlem 

3 . Fluid Dynamics 

4. Test Data Analysis 

5 . Summary 

6 . References 
Appendix A 
i^endix B 


page 

i 

iii 

iv 
1 
k 

44 

53 

60 

64 

65 
73 




IHustrations 


Figure 

1 

Mathematical Model ajid Coordinate Systems 

6 

Figure 

2 

Exact and Approximate Nutation Angle Time 
Histories for the Nutation Synchronous Jtode 

12 

Figure 

3 

vs e 
cn 

13 

Figure 

4 

T vs T1 
cn ' 

l4 

Figure 

5 

cn 

15 

Figure 

6 

T vs ff 
cn 

16 

Figure 

7 

T vs 7 
cn ' 

17 

Figure 

8 

T vs 6 
cs 

20 

Figure 

9 

T vs n 

cs 

21 

Figure 

10 

T vs h 
cs 

22 

Figure 

11 

T vs 0 
cs 

23 

Figure 

12 

Tc3 vs 7 

24 

Figure 

13 

T and T vs g 
cs cn 


Figure 

l4 

T vs Ou 
cs 2 

39 

Figure 

15 

T vs 01 
cn 2 

43 

Figure 

l6 

Normalized Energy Dissipation Per Cycle 
as a Function of ^ 

52 


- ill - 



A, B, 


C 

H 



R 


a 

h 

°d 

g 

g 

h 

k 

m 

t 

6 

5 

€ 

7 

7 

T) 

V 

<■ 

V 

0 

SB 

“t 


Sj-ribols 

- transverse moments of inertia of the satellite 

- spin axis moment of inertia of the satellite 

- angular momentum vector with con5)onents 

- transverse angular momentum 

- moments and products of inertia of the fluid slug 

- radius of the annulus 

- radius of tube 

- h/R 

- viscous damping coefficient 

- gravity 

- g/Rii^ 

- height of annulus above satellite center of mass 

- (sin (r/2))/(7/2) 

- mass of fluid 

- dimensionless angular velocity coniponents 

- time 

- offset of center of annulus from spin axis 

- 6 /R 

- xaS^/ky 

- angle of fill of fluid in annulus 

- fraction fill, y = y/2j( 

- diioenslonless dashing coefficient 

- angular position of offset in transverse plane 

- kinematic coefficient of viscosity 

- initial spin rate 

- angular velocity with con^onents o^, 

- transverse angular velocity 
dimensionless time, r « flt 
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0 . nutation angle 

4 - precession angle 

\|r - Euler angle 

B - angular position of slug 
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1 . Introduction 


A ring partially filled with a viscous fluid such as mercury 
was one of the first nutation dampers used on spinning satellites. 

The first analysis of the partially filled viscous ring danqper was 
performed hy Carrier and Miles * . They assumed that the motion of 
the daB 5 )er did not appreciably affect the precession rate of the sat- 
ellite hut acted only as a source of energy dissipation. With this 
assumption the motion of the fluid in the tube was then treat i. as 
a fluid mechanics problem and an approximate solution to the Navier 
Stokes equations was obtained. Their solution showed that the fluid 
behaved as a rigid slug for a nutation angle greater than one degree. 
At these large nutation angles the problem was then treated as bound- 
ary layer flow over a flat plate with the width of the plate being 
equal to the perimeter of the tube. However this analysis did not 
completely treat the problem as there are two distinctive modes of 
motion for a nutation angle 9 greater than one degree, and in one 
of these the fluid does not behave like boundary layer flow. The 

■a 

next analysis was performed by Carwright ' , et. al., in which they 
assumed the fluid mass behaved like a particle of equal mass moving 
in a tube with a viscous damping force. Their analysis revealed that 
there are two distinctive modes of motion which they called the nuta- 
tion synchronous mode and the spin synchronous mode. Although there 
were some minor errors in their equations of motion they correctly 
analyzed the nutation synchronous mode but failed to analyze the spin 
synchronous mode. 

Interest was revived in this problem when the failure of the 
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ATS-5 satellite was attributed to the energy dissipation caused by 
fluid motion in the heat pipes . Consequently it was desirable to be 
able to predict more accurately the energy dissipation in this type of 
damper. Also other satellites such as Helios will employ partially 
filled rings for nutation dampers but during a portion of the flight 
the satellite will be spinning about an axis of minimum moment of 
inertia. Since this is an -unstable configuration it is ’"ery important 
■to be able to predict the rate of energy dissipaton. 

5 

As a resiilt Alfriend approached the problem in the same manner 
3 4 

as Cartwright ^ and obtained equations which approximately describe 
the motion in both the nutation and spin synchronous modes . The prob- 
lem -with this approach is that it must be accompanied with a method for 
calculating the dauping constant. Leibold^ suggested assuming steady 
flow in straight pipe as a means of calculating this danping constant. 
Howe-ver there is an error in his equation describing the motion of the 
satellite . 

In this s-tudy the fluid is assumed to behave as a rigid slug but 
but since the fluid may fill up -to 50^ of -the ring it is asstaned the 
fluid is a rigid slug of finite length, not a particle. 

In Section 2 a description of the mathematical model is given and 
for a symmetric satellite approximate solutions are derived for the nuta- 
tion angle time history and corresponding time constants in both the 
nutation synchronous and spin synchronous modes. These approxima-te 
solutions are then compared to those obtained by numerical integra- 
tion of the exact equations of motion. In Section 2,2 the effect of 
a ffTTw-ii offset of the center of the ring frcm the spin axis is inves- 
tigated. This is necessary since an offset of l/4” is planned on -the 
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Helios satellite. An analysis of the test data was a part of this 
investigation so the effect of gravity on the nutational behavior is 
Investigated in Section 2.3. An investigation of the effect of a stop 
in the tnbe is investigated in Section 2.4, In Section 2.5 the results 
of Section 2,1 are extended to the asynmetric satellite. 

Several possible methods for determining the damping constant 
are given in Section 3. An analysis of the test results is presented 
in Section 4. Finally a summary of the results and conclusions are 
given in Section 5 




2. Statoment of the Prohlem 


In the presentation of the analytical results the simplest prohlem, 
vhich is the symmetric satellite with no gravity and no ring offset, is 
solved first. The effects of gravity and ring offset on the dymmetric 
satellite prohlem are then detemined. Finally the solution of the 
asymmetric satellite with no ring offset and no gravity is given. The 
symmetric prohlem is considered first rather than solving the asym- 
metric prohlem and sic5)lifying the results for the symmetric case be- 
cause hy solving the symmetric prohlem first one gains more insight 
into the prohlem. 

The mathematical model is an asymmetric rigid body (satellite) with 
principal moments of inertia A, B and C and corresponding principal 
axes represented hy the x, y and z axes shown in Figure 1. The z 
axis is the spin axis. A tube of radius R is attached to the rigid 
body at the point (hcosv, 8sinv, h>. 6 is the offset of the center of 
the ring (tube) from the spin axis. Moving in the tube is a rigid 
slug of mass m which fills a portion of the tube, the angle of fill 
being 7. The other assuii5>tions in the development of the equations 
of motion are l) the center of mass of the system and the center of 
mass of the satellite are coincident, 2) the friction force on the 
fluid slug can be represented by a linear viscous force, and 3) gravity 
acts only on the fluid slug. The first assuniition is made because it 
simi>lifies the equations of motion considerably and the effect of the 

O 

motion of the system center of mass is negligible since it is of 0(e) 
where e is a small parameter which is defined later. Also, in tests 
the satellite center of mass is a fixed point. The third assumption is 
made since gravity has an effect only in the tests and the satellite 
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without the fluid was statically balanced before the tests . 

The equations of motion which are derived in Appendix A are 


r (1+a ) (1-a ) 
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The coordinate systems and angles are defined in Figure 1, 
p, q and r are the coirponents of the dimensionless angular velocity 
along the u, v and z axes and 0 defines the position of the slug 
in the tube. The and ' are the moments of inertia of the sluiz and 
their derivatives are given in Appendix B. The independent dimensionless 
parameters of the system are Og, e, h, n, 7 or y, g, 5 and v 
where 

= C/A 

cTg = C/B 

o ^ 

e => mR /A7 

h = li/R (2.6) 

Ti * c^/mfl 

g = g/Rn^ 

S » 5/R 

Primes denote differentiation with respect to the dimensionless time 
T ■ At where A is the initial spin rate, 0, t and ^ are the Euler 
angles with 9 being the nutation angle and ^ the precession angle. 

The Euler angles are detezmined from 


0» 

• p cos(tt^) - q sin(iH^) 

(2.7) 


» [p sin(irt^) + q cos(\|f4^)]/sin0 

(2.8) 


• r - 4’cos0 

( 2 . 9 ) 
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2.1 Symmetric SateiJLite with Zero Offset and Zero Gravity 

Since there are no external forces the angular nomenttm is con- 
stant. Letting the reference direction he H (the vertical when 

~ z 

gravity is present) the nutation angle 9 can be determined frcan 


h 

tan9 = ^ 
xl 

z 


( 2 . 10 ) 


instead of using (2.T). is the spin axis component of the angular 
momentum and is the transverse component, i.e.. 


. 

t u V 


( 2 . 11 ) 


Substitating for and and expanding in a power series in 


e gives 


tan0 = — + 0 (e ) 
c ' 


( 2 . 12 ) 


where 


2 2^2 
03^ = p + q 


cr = cr^ = cr^ 


(2.13) 


However to use (2 .10 ) to obtain the variation of 9 one would have to 
obtain p, q and r through 0(e)^ which is no easy task. Rather 
than using (2.10) or (2,12) one can obtain a good approximation of 
9 by differentiating (2.10) and then integrating the restilting 
differential equation. Differentiation of (2.10) gives 


H‘ H' 

0» = i 

H H. 

z t 


pH -qH 
^ V ^ u 

H, 


(2.1U) 


:r 



- 9 - 


Substituting for and and expanding in a power series in e yields 

0' » [^p+(l4^’)t> sin(r/2)]^+0(e^) (2.15) 

The advantage of (2.15) is that to deteimine 0 to 0(e) one only 
needs the first approximation of p, q and p . 

Damper Motion 

A symmetric rigid body which is spinning about its axis of symmetry 
has a constant nutation angle when no damping is present. The transverse 
angular velocity vector rotates at a rate of cDcos0 and the body ro- 
tates relative to at a rate of (l-o)n. When no daii 5 )ing is present 
the center of mass of the fluid slug will be flung outward as far as possi- 
ble which will be along cn. or the plane formed by H and the z axis, here- 
after called the nutation plane. The fluid slug will then be moving at a 
constant rate of (l-cr)n with respect to the body. Introduction of a small 
amount of damping causes the center of mass of the fluid slug to move off 
the nutation plane to an equilibrium position where a component of the cen- 
trifugal force balances the friction force. This type of motion is called 
"nutation synchronous” motion . In this mode the fluid slug is moving at 
a constant rate with respect to the body, hence the energy dissipation rate 
is a constant. If o > 1 the nutation angle decreases which causes a decrease 
in the centrifugal force and the fluid slug center of mass moves further from 
the nutation plane. Eventually the centri‘**ugal force is not large enough to 
balance the damping force and the fluid slug begins to be dragged around with 
the body while oscillating in the tube . This type of motion is called spin- 
synchronous motion . 

The purpose now is to determine the behavior of the nutation angle in 
these tvTo modes as a function of the dimensionless parameters e,Ti,b, o and y. 
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Nutation Synchronous Mode 

Ejqjanding the equations of motion, Equations (2.1)- (2 .4), 
power series in e and dropping all terns of 0(e) gives 

in a 

II 

U 

(2.16a) 

P' + (J^-^’)q = 0 

(2.16b) 

q* - (>^-p')p = 0 

(2 .16c ) 

P” + t) 3» + -^^q + boq(sin(7/2))/(7/2) = 0 

(2.l6d) 

where ^ = 0 - 1 . 

Letting 

a = p - Xt 

the solution to (2.l6b) and (2.l6c) is 

(2.17) 

p = -03^ cos a 


q = 03. since 

(2.18) 


where it has "been assiamed that q = 0 and p = 0 at t = 0, Thus a measures 
the _position of the center of mass of the fluid slug with respect to the 
nutation plane. Substituting (2.1?) and (2 .18) into (2.l6d) gives 

CD. cos (7/2 )cos a 

a” + Tja’ + hocD^ sina (sin(7/2))/(7/2)[l (2.19) 


A particular solution of this equation is 

tanScos (7/2 )cos a 

.inajl £-] 


a = a where 
e 

= 

ha^ tan0sin(7/2 ) 


( 2 . 20 ) 


whesv* 03^ = ctanS has been used. Thus the fluid slug will remain in 
the nutation synchronous mode as long as (2 .20) is satisfied. Once 0 
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■becomes sTtiH-li enou^ so that o = ± it/2 the fluid slug goes into the 
spin synchronous mode. The transition angle 0^ from one mode to 
t)ie other is 


_ Ti|g-l| 

t " 2 sin(7/2) 

Substituting for p and q in (2 .15 ) gives 


( 2 . 21 ) 




[ - (TtanQcos (y/2)cosa 

^ h(<W') 


~sinasin(7/2) + 0(€^) (2.22) 

« 


Substituting a = a and dropping terms of ©(e'^) one obtains 

© 


tan0 0» * - iDi - C — j j. 

CT 


(2.23) 


for which the solution is 


where 


COS0 = COS0 exp(r/T ) 
o ' cn 


g 

eri7(g-l) 


(2.24) 


(2.25) 


Thus in the nutation synchronous mode the cosine of the nutation angle 
exhibits exponential behavior. If there are N daB5>ers the time con- 
stant is 



(2.26) 


It has not been assumed that 0 is small thus (2.24). is valid for 
0 < 0 < jt/2. For small 0 the nutation angle time history is 
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e 



2t 


r ^ 
cn o 



(2.27) 


A comparison of the approximate solution given hy (2.24), to the 
exact solution is shown in Figure 2. A conparison of the time constant 
given hy Equation (2.25) and an exact time constant ohtj.ined hy inte- 
grating the equations of motion, Equations (2.1)- (2 .4), is given in 
Figures 3-7. 


Spin Synchronous Motion 

Substituting (2.18) into (2*l6d) yields 


ovff 0) cos(7/2)cos(S-Xt) 

+ ^a3^sin(7/2)sin(p-XT)[l--i ^ ] 


( 2 . 28 ) 


Assuming that « 1 and using the first iterate of a Picard iteration 
scheme a good first approximation of the steady state solution of Equa- 
tion (2.28) is 


e 


ho^ 


tan 0sin(7/2)/(r/2)[E sin(S^-^r)+F cos(S^-^t)] (2.29) 


where" 




F = 



The nutation euigle differential equation is 


i 0* « cos(s->^T)+(l-»fl' )l3 sin(7/2)]sin(p-3w-) . (2.30) 

K 7 ^ 

I 


I The exact time constant is obtained hy integrating the exact equations 

over a suitable period of time, assuming e:q)onential behavior for cos6 
and calculating the time constant. 


#■ 
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p 

Assuming that 9 is small enough so that terms of 0(0 ) can he neglecte i 
and assuming that the change in P is small so that sinG0 = bp and 
cos 6 P =1 ( 2 . 30 ) reduces to 


where 


0* + 0( + K cos2^t+k sin2^T ) 

1 2 
cs 




- ~ sin O' /2 ) sin^T ( 2 . 31) 

JT 


( 2 . 32 ) 


The solution of ( 2 . 3 I) is an infinite series but the and terms 
contribute nothing to the exponential decay of the solution. Thus 
the lB 5 )ortant part of the solution is 

9 . 9 sinXT+XT® (coaV-e''"'’’”)] ( 2 . 33 ) 

° +1) “ “ 

C5 

Comparison of t given by (2.32) and the time constant obtained from 
c s 

numerical integration of the exact equations of motion is shown in Figs .8-12. 
2.2 Damper Offset Effect 

The purpose of having the center of the ring offset from the spin 

aarf-s is to guarantee that the fluid will act as a rigid slug. The 

1 2 

analysis of Carrier and Miles * showed that for very small nutation 
angles the fluid would spread out along the outer wall of the tube . 

However, if there is an offset the centrifugal force will be greatest 
(for very gwaJi nutation angles) in the direction of the offset and 
the fluid will '.end to lump there. Consideration of the physical 
situation is an aid in determining the effect of the offset. Let 
the offset angle be the angle between the spin axis and the vector 
frcm the satellite center of mass to the center of the ring. For 
r na ll offsets the offset angle is smaller than the 
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transition angle. When 6 > 9,j, and no offset the fluid slug maintainr 
a fixed position with respect to the nutation plane since the centri- 
fugal force is constant. The effect of the offset should be to cause 
small oscillations of the slug about this equilibrium position of the 
slug. The result would a small change in the rate of energy dissipa- 
tion which would cause a small change in t . ■ Numerical integration 

cn 

of the equations of motion has verified these conoectures . For the 
Helios satellite the offset of ifh" caused less than a 10^ change in 
^cn* ^^cn since the rate of energy dissepation increases.) 

In the spin synchronous mode (0<0^) and zero offset the fluid 
slug moves slowly around the tube while oscillating. Since the offset 
would create a point in the ring where the centrifugal force is a 
maximum the fluid slug should oscillate about this point instead of 
moving slowly around the tube with the result that the change in t 

cs 

should be minimal. Again numerical integration of the equations of 
motion verified these conjectures. Except for one cane the change in t 

CS 

was less than 10^ for the Helios satellite with a l/4" offset. The one 

exceptional case will be discussed in Section U. 

2.^ Gravity Effect 

Since the satellite was bal anced without the fluid it is assumed 
that the gravitational force acts only on the fluid slug. Following 
the procedure of Section 2.1 the approximate equations of motion are: 


r = 1 (2.34a) 

p* + (J^-p‘)q = 0 (2.34b) 

q«-(X-g*)p = 0 (2.34c) 

P'* + T)p‘ + + bokq + kg sin0cos(P + t) ® 0 (2.34d) 

7 
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where 6 is the angle between the spin axis and local verti'.al. Note 
that the only term involving gravity is the last term of (2.34d). The 
solution of (2.34b) and (2.34c) is 

p = -03^sin(V-p) 

:2.j5) 

q = -o^cos (^T-3 ) 


Since gravity is present the angular momentum is no longer con- 
stant, consequently H does not coincide with the vertical. Let 
d be the angle between H and the spin axis . Then 

A ^ 


tan 0, 

A 



(2.36) 


Differentiation yields 

®h * €7C^^^P+hk(l+S‘ )]sin(S-XT) -tcyk g sin(p-?^T)cos0j^cos (0-0^^) (2. 37) 

where 


i|r = -Xt - rt/2 


(2.38) 


has been used. 

The assumption is now made that the change in is small ccmpared 

to the change in 0, and that 0 = 0^^ at t = 0. This allows one to use ® 
in the equations . Numerical integration of the equations of motion shows 
that this is a reasonable assumption. Equation (2.37) becomes 


01 _ e 7 [^^Ip+bk(l 4 ^’)+kg cos 0 ]sin(P-Xr) (2.39) 

Nutation Synchronous Mode 
As In Section 2.1 define 

a = p - Xt (2.4o) 




Substitution into (2«34d) gives 



a" + 


Tia'+ [cD.lDck+gk sin0 




* 


(2 .41) 


There is an eqaillhritm value a = which is given hy 

[h(j\ tan0(l + -2^os0) - tan^ e ^ - ^^^ cos a ]sina = -ill (2.42) 

b(T^ 1 ® ® 

The transition angle between the two modes is obtained by setting 
= ± jt/2, which gives 


tan0„(l + -^cos e ) = 

b(i bcTc 


which for small 0 becomes 


aisil 

bAo 


where 

G = 1 + g/(T^ 

Siibstituting (2.42) into (2.39) yields 

tan0 0« = iMllii 
a 


(2.43) 


(2.44) 


(2.45) 


(2.46) 


which is the same equation one obtains without gravity. Thus the 
only effect gravity has in the nutation synchronous mode is to change 
the transition angle 0^. The solution to (2.46; is 


COS0 = COS0 e2q>(T/T ) 
O Cu 


(2.47) 


where 


£711(0-1) 

Spin Synchronous Mode 

Using the procedure of Section 2.1 an approximate solution of 
(2.34d) for small 0 is 

&P = K(lsin(p^-lT)-iicos(p^-lT)) 


(2.48) 


(2.49) 
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where 


K = 


kho^Qg 

X(XV) 


(2.50) 


Substituting into (2.39) ancL using sin5P= 5S and cosSp= 1 gives 


where 


0* + ^ — 0 = -€ 7 hk G sinXr 

cs 


r - 2(g-l)[(g-l)^+Ti^1 


(2.51) 


(2.52) 


Thus gravity can have a substantial effect on the nutation angle in the 
spin synchronous mode as the time constant is reduced by a factor 1/G . 

The approximate time constants given by (2,48) and (2.52) are com- 
pared with those obtained by numerical integration in Figure 13 . It is 
seen that the agreement is very good, 

2.4 Effect of Closed Ends in Tubes 

The manufacture of heat pipes is siii 5 >ler if there is a stop (closed 
end) in the pipe. This analysis has been luidertaken to determine if it 
is advantageous from the standpoint of nutation damping to have a stop 
in the tiibe. The physical situation is so ccn 5 )lex that it defies an 
accurate analysis , What happens to the fluid when it impacts the stop 
is not known, thus a very simplified analysis has been performed to 
try and determine if the stop increases or decreases energy dissipation. 

It is assumed that the fluid slug is a point mass and that its motion 
is governed by (2.l6d). Since the nutation angle is essentially constant 
over one cycle it is also assumed that the nutation angle is constant. 
Equation (2.l6d) is then integrated over a number of cycles taking into 
accou::t c'llisions and the average amount of energy dissipation is then 
determined. This is then compared to the change in kinetic energy for 
no stops which can be computed analytically. 
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For a continoiis tube the rate of energy dissipation is 

• 2^2 - 2 
T = P = 


(2.53) 


The change in kinetic energy over one cycle is 
= 


T 

/t dt 


where is the period and 




, g» , 

|a-i|n 


This gives 


- o2 






where P* is the average value of P’ . Let 


.2 


2^ = 


Zsrp 


A£7fl 


then 


In the nutation synchronous mode 

P* = (1-<t) 


(2.54) 

(2.55) 

(2.56) 

(2.57) 

(2.58) 


hence 


2^ = 2rtTija-l( (2.59) 

In the spin synchronous mode. 

3' ' [V:os(g ->-T>fT|eto(p ->t)] 

(^V) ° ° 


hence 


2St * 




(2.60) 


I a-l| (J^W) 

The change in kinetic energy when there is a stop in the pipe is 
obtained by numerical integration of (2.l6d) and 
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• - 

^ = •np'^ (2.61) 

To this ve have to add the loss in kinetic energy from the collisions 
which is determined by assuming the slug collides with an object of 
Infinite mass . The loss in kinetic energy from a collision is 

c 2 "^i 

where e is the coefficient of restitution and is the angular 
rate just prior to collision. 

Runs were made for both a stable and unstable configuration in 
the nutation synchronous mode and spin synchronous mode. Table 2.1 
gives the results for several values of the coefficient of restitution e. 
As one can see the change in kinetic energy is relatively independent 
of e. The increase in energy dissipation is significant in the nuta- 
tion synchronous mode but small in the spin synchronous mode. This 
is reasonable since in the nutation synchronous mode the angular rate 
of the slug with respect to the tube is |a-l| whereas in the spin 
synchronous mode it is small. In tests run by Hu^s Aircraft^ ^ on 
the AIS-V heat pipes there was a increase in the change in kine- 
tic energy in the heat pipes with stops . No con5>arison can be made 
with the Hughes results because the damping constant tj is not known. 
However one can conclude that putting stops in the rings can cause a 
substantial increase in the energy dissipation which will result in a 


decrease in the time constant. 



Table 2.1 


In all runs = 0.1 


Nutation Synchronous 

ff = 1.2 , 0 = 6" , b = 

3 

e 

^(no stop) 

AE( 1 stop) 

^ increase in K.E. 

0.2 

0.126 

0.969 

670 

0.9 

0.126 

1.002 

700 

0.6 

0.126 

1.045 

730 

0.8 

0.126 

1.060 

74o 

1.0 

0.126 

0.839 

565 


Nutation Synchronous 

iT * 0.8 , 0 = 6" , b = 

3 

e 

^(no stop) 

^(1 stop) 

^ increase in K.E. 

0.2 

0.126 

0.377 

200 

0.4 

0.126 

0.391 

210 

0.6 

0.126 

0.403 

220 

0.8 

0.126 

0.431 

24o 

1.0 

0.126 

0.520 

313 


Spin Synchronous a ® 

1.2 , 0 = 1" , b = 0.8 


e 

At (no stop) 

^(1 stop) 

io increase in K.E. 

0.2 

0.0127 

0.0162 

27 

0.4 

0.0127 

0.0163 

28 

0.6 

0.0127 

0.0166 

31 

0.8 

0.0127 

0.0163 

28 

1.0 

0.0127 

0.0178 

4o 


Spin Synchronous a ~ 

0.8 , 0 = 1' , b = 0.8 


e 

^(no stop) 

^(1 stop) 

increase in K.E. 

0.2 

0.0039 

; 

0.0048 

23 

0.4 

0.0039 

0.0045 

15 

0.6 

0.0039 

0.0045 

15 

0.8 

0.0039 

0.0048 

23 

1.0 

0.0039 

o.oo4o 

3 









2.5 Asynmetric Satellite with Zero Offset and Zero Gravity 


In developing the approximate solution for the asymmetric satellite 
it is advantageous to tise the components of the angular velocity along 
the X and y axes rather than the u and v axes . Iiet p and q he the dimen- 

^ y 

slonless components of the angtilar velocity along the x and y axes, then 


p = p cosp - q sinp 

3C 


(2.62) 




= p sinP + q cosP 


Using the same procedure as in Appendix A the equations of motion become 

Px-*-^V°12V ' "7 ( — 2 "^V p^tq^r) 


(1 -i ) 

+ [ (-p^-qyr.2q^ ' )cos2p 


+(-q^+P/+2p^’)sin2p] - 

+ (r *4^")cosP-(r4p ’ )^sinP+qy (p^cosP+q^sinPlI 


(2.63a) 




(1 -I ) 

+ - ^g - — [ (+q^-P^r-2P *p^)co82P 


+(-Pi -qyr-2P’qy)sin2p] (2.63b) 

+ i P^(r^^’) + I r(r'-»^")8inp+(r4P')^co8P 

ZZr X UZ 
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air‘ + (a^-l)qyP^ = iyyp' 


( 2 . 63 c) 


P" + nP’ + — _q^^+ 

+ ' [(p^-q^)sin2p-2p^q^cos2p] = 0 


(2.63d) 


The method used to obtain the* behavior of the nutation angle in 
the symmetric case can still be applied in the asymmetric case but the 
resulting expression for 0*, Equation (2.l4), does not sinplify to 
0' = 0(c) because in the asymmetric case 0 is not constant when there 
is no damper. The result is that to use (2 .14) higher order approxi- 
mations of ^a^d qy would have to be developed. To get around this 
difficulty a variation of parameters approach will be used. When there 
is no darker the nutation angle 0 oscillates^ thus the time constant or 
damping constant obtained will measure the increabe or decrease of 
the maximum value of 0 at the end of eaoh oscillation. 

For e = 0 the equations of motion for the satellite are 


fi * <v^i2 V ■ ° 




^ qj 
o^erg Va 


( 2 . 64 ) 


The solutions for p , q. and r an elliptic functions . It iriU nesr be 

* y 

assumed that the natation angle is small enough or the asynoetry is 
enough so that the elliptic functions can be r^tlaced by tbe 
first term in their trigonometric esqpansions. This will limit the 
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approxijiiate ?;olutlons to small angles. The approximate solutions 
when (or^-l) and have the same sign are 

Vy. = -oi^cos (J^T+x) 

Qy = -oi^M sin(V+x) (2.65) 


r = 1 

where x is a phase angle end 



X® « (0^-1) (Og-1) 


The n. tation angle for small angles is given hy 


. . 

H 

z 


f .2 2^_2 2 ^ 

(A 05 +B 05 ) 

— Z 

^CD 


1/2 


+ 0(e) 


In tezms of the dimensionless variables this becomes 


, 2 2 2,^/® 


Substituting (2.65) for p^, and r gives 


1(L^+o.-2| Iv®iI 

9 - [—=-5= + — ■=5-^cos2(>^t+X)] 


l/2_ 






77^ 


(2.66) 


(2.67) 
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A variation of parameters approach will now he used to ohtain o)^ 
which rfill then give the approximate nutation angle time history. 
Assuming and •)( to he functions of time, substituting for 

Px’ V 

linearizing with respect to and solving for o>^ gives the following 
differential equation for co^. [We do not need to solve for;;^'] 


and r* into (2.63a,h) neglecting terms of 0(e ), 


tn; 






+ 2 


*^1 “ 


(i ~f ) 

~ gg" ^1*^2 ^ ^ ^l'*' °2~^ ^1 "^l" °2 ^ 


(I -i ) 

+ - ^2 ' ff ^ ^ ^ '^l ° 2 ~ °l" °2 ^ V “^l ) 


( 2 . 68 ) 


(i -i ) 


^ ^ -(a„-o, )P'cos2(J^t+X) 


2or^ '2 1' 


+ (p”cosP- (l-t^ ’ )^sin3 )co5 (lT+x)+(P”sinp+(l4^ ’ )^cosP ) 


+ (l4§*)^J^sin2(XT+X)/ 


Spin Synchronous Mode 

Substituting (2.65) for p^ and q^ into (2.63d) and neglecting 
2 

terms of 0(o3^) gives 
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P" + Tp* = 


UZ 


• 

(a -1\ 



iV^I 


\ 1 


cospsin(^T+X) 


(^T+X) 


(2.69) 


CO, 


As in Section (2.1) an approximate steady state solution is 

developed "by using the first iterate of the Picard iteration vhich 

is just the solution of (2.69) with P = P on the right-hand side of 

o 

(2.69). 

gp = p - P = ILcB.cosp. (Ticoslr+Xsin^r) + K„a).sinp (+tisin^T->coslT) (2,70) 

O 1 t O 2 t O 

a/2 


K, = - 


Kg = - 




„^X(xV) 


(2.71) 


ihus the first approximation of p gives an oscillation about p^ 
with the magnitude of the oscillation proportional to cu^. 

It has heen assumed that cb. is small, hence we can use 

sinP = sinp + gpcosP 
o o 


cosP = cosP - gpsinP 
o o 

2 

Using these approx3mations and neglecting terms of 0(o>^) after sub- 
stituting for P and p* in (2.68) gives 
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co^ + :~+ K^sin2(J^T+;^)+»t2COs2(V+;;^) oi^. = eyCic sin(XT+;^)+ic^siii2(^T+;^)] 

cs ^ 

(2.72) 

^ere the are constants and 


_iz 


7b\^T 


cs 


2 a^(Ati^) 


f"lVV°l> . 2. . 




(2.73) 


As in Section (2.1) the solution of is an infinite series hut 
the important part of the solution is the exponential decay given hy 




- -t/t 
C0j.e ' cs 


(2.74) 


The solution for the nutation angel 6 hecomes 


e = 


|0i+0‘2-2| Iv°1^ 


'“o -T 1/p 

+ — ^ -3- ^ - cos 2 (Xt+x ) ^ -r e 


cs 


(2.75) 


Thus the maximum value of 0 in each oscillation decays with a time 
constant of t given hy (2.73)* Note that t is a function of the 

OS OS 

position of the slug in the tube as shown hy the prescence of in 

(2,73)* N^^merical integration and the second iterate in the Picard 

iteration solution of (2.69) reveal that the slug moves slowly in the 

tube in addition to its oscillation. The result is that t is a 

cs 

slowly varying function of time. For a design criteria one should 
use the value of which gives the maximum time constant, 

Conparison of (2,73) with an "exact" time constant obtained from 
the exact equations of motion is given in Figure l4. The conparison 
is made for P^ = 0, In the numerical integration a very small offset 
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was included in order to make the slug oscillate about p=0. The "exact" 
time constant was obtained by assuming exponential behavior for the 
maximum value of 0. From Figure l4 we see that the time constant given 
by (2.73) is a good approximation of the exact time constant even when 
(ffg-cr^) is not small. 

Nutation Synchronous Mode 

■When the slug maintains a fixed position with respect to 

the nutation plane in the nutation synchronous mode. As the nutation 
angle slowly changes the position of the slug changes slowly. In the 
asymmetric case we would expect the behavior to be similar but in the 
asymmetric case 9 oscillates . Thus for small oscillations in 9 we 
would expect the slug to oscillate about some equilibrium point. For small 
oscillations of the slug the energy dissipation resulting from this 
oscillation would be small with respect to the total energy dissipa- 
tion. The result should be a solution like the one for the symmetric 
case. However as increases or for large values of 9 the magni- 

tude of the oscillations in 9 increases with the result that the energy 
dissipation from the oscillations of the slug could became significant. 

In this case we wo'uld expect a motion which is a combination 

of the nutation synchronous and spin synchronous modes of the symmetric 

case. 

As in the symmetric case let a define the position of the slug 
with respect to the nutation plane. 

p * V + X + a (2.76) 


Substituting in (2.69) gives 
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^ “t * 

a" + TP' + t(a^ff2+(T2-ai)M +(a^cT2+o-^-a2)]sina 


+ [(a^ff2+a^-a2)-(o'^q2+ar2-a^)M*]sin(2XT+2)(;+a ) = -tjX 


* a -1 1/2 
where M = (- - t) 

2 


The coefficient of sin(2lx + 2X +ot) is small conrpaxed to the coeffici- 
ent of sin a, in fact it vanishes as for ^ approximate solu- 

S"t 

tion for a obtained hy the 1 iterate of a Picard iteration is 


a = a + 
e 




[2lsin(2lT+£?(;+a )+ cos ( 2 lT+ 29 (+a )] 
€ 6 


where 


sin a = 


-2<j^nl 


e " ai^[ jM*+(a^o- 2 +a^-a 2 j] 


(2.79) 


Thus the magnitude of the oscillation is proportional to o)^ or 0 and 
|ffg-(T^|. This oscillation in a will he neglected in deteimining 
the solution for cu^. Substituting a = in (2.68) yields 

1 (1+1)2 ^ 

cuJ = e/ (l-rtl )sina + oscillatory terms (2.8o) 

^ 2M ® 


The oscillatory terms contribute nothing to the decay of so they 
will be dropped. Substituting for sin gives 


Vt “ “ -ern(vi) 

M 


(2.78) 
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which has the solution 


= [oj^-2e7ii(T^(ar2-l)T]^/^ 


(2.81) 


Substituting into (2.6?) the solution for 0 becomes 




1/2 [a?- 2 € 7 TlcT (a -1 )t]^/^ 

2 1/2 ^2.82) 

[2ct^(02-1)] ^ 


Letting 9^ be the initial value of 0 and ‘assuming it is the maxi- 
mum for the first oscillation we get 


0 = 


I a^+Og-a 1+ { Og-CT^ I cos2 (J^T+x) 

1 W^l.^^°2’°'ll 


1/2 


(02.2T/T 

o ' cn 


(2.83) 


where 


cn 


2(j^sign((j^+(T2-2) 
e^Tido-^+Og-al + lag-a^l) 


(2.84) 


A comparison of the dancing constant given by (2.84) and an "exact" 
obtained by numerically integrating the exact equations of motion 
and assuming the maximum value of 0 during each oscillation behaves 
according to 


0 

max 



is given in Figure 15. Figure 15 shows that the approximate solution developed 
is a good approximation even when (<^2**°i^ small. 




3. Fluid Dynamics 


In Section 2 approximate equations were developed for the nutation 
angle time history and the corresponding time constants . Comparison of 
these approximations with exact solutions obtained from numerical in- 
tegration showed excellent agreement. However, this means that the com- 
parison is good for the mathematical model of the system. Two import- 
ant questions which still need to be answered are: l) How good is the 

mathematical model? 2) If the mathematical model is valid how does one 
calculate the damf'ing constant? The answer to both of these questions 
can come only from testing. In this section several methods for calcu- 
lating the daurping constant tj for symmetric satellites are presented. 

The analysis of some test res\ilts is given in Section U. 

Nutation Synchronous Mode 

In previous studies two approaches have been used to calculate the 
damping constant. Both approaches have drawbacks in that certain assymp- 
tions are made which are not completely valid. One approach is to aiodel 
the motion of the fluid as steady flow in a pipe, and the other approach 
is to model it as boundary layer flow over a flat plate with t' e width 
of the plate being the perimeter of the pipe. Both of these approaches 
must be considered for both laminar and turbulent flow. 

a) Steady Flow in a P_pe 
1 . Laminar Flow 

The development of steady flow in a strai^t pipe can be found in 
sQmost any standard fluid mechanics text such as Reference 7 . 

Solution of the Navier-Stokes equations with a flux of 


.kk. 
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(i = na 


(3.1) 


where R3 is the velocity of the fluid slug relative to the ring and 
a is the radius of the ring, gives 


u = 2B^(l-(r/a)^) 


( 3 . 2 ) 


The shear stress at any point is 


5u 

Where p, is the viscosity. 
Thus the total viscous force is 


( 3 . 3 ) 


F = (2ita)(R7)T|^^ = 




Equating this to the force detenained from the dynamic analysis, i.e.. 


F * cv = cK^ = Tj m li R p 


( 3 . 5 ) 


gives 


Tl = 8(-f-) 


( 3 . 6 ) 


where v is the kinematic viscosity. 


The quantity (^ — ) is a Reynolds number hut is not the standard Reynolds 


number for this type of analysis . The standard Reynolds number is 


«e = 


( 3 . 7 ) 


Equation (3.6) holds for R < 2000 since the flow is lamintr for R <2000. 
The analysis above assumed flow in a strai^t pipe hut we have flow 

in a curved pipe . The correction factor for flow in a curved pipe as 

8 

given hy Schlicting is 
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t/t^ = 0.1064[R^(a/R)^/^]^/^ (3.8) 

where is the shear stress in the straight pipe and t is the shear 
stress in the curved pipe. Equation (3.8) is valid for 10^’^ < (a/R)^'^^ <10^. 
For Rg = 2000 and (a/R) = l/lOO the increase in shear stress is 50^* 

Therefore one should take into account the curvature of the pipe. Thus 
for laminar flow the dan^jing constant 'becomes 

n = (3.9) 

a n 

The as sung) ti on here is that we have steady flow in a pipe but a 
certain length is required for steady flow to develop. For flow from 
a cistern into a pipe this length (Ref. 8 pg. 301) is 

6 = 0.0575 a R^ (3.10) 

which for R^ = 1500 is 5 = 86a. But "the length of the fluid in many 
cases may not be much more than 86a. Thus the length of fluid required 
for steady flow to develop may be about equal -to the lejagth of the fluid. 

Thus -there is an error in assuming steady flow. 

2 . Turbulent Flow 

Blasius (Ref, 8 pg. 339) developed for -the shear stress for stea.dy 
turbulent flow in a straight pipe the enpirical result 

T - 0.0791 R'g'^(|p 4) (3.11) 

where u is the mean velocity R^ . This result is valid for <10^. 
m e 

Using the same procedure as before -to calcula-te the dancing constant 



one obtains 


- 47 - 


1 

a fl ^ 


(3.12) 


The correction factor to take into accotoit the fact that the flow is 
in a curved pipe is 


t/t^ = 1. + O.C75 


(3.13) 


The correction factor for turbulent flow is smaller than that for laminar 
flow and can be neglected since it is usually less than 10^. With the 
correction factor the damping constant becomes 



If the flu'^.d is free from disturbances at entry the flow in a smooth 
pipe for some distance 6 from the entry will be laminar even though 
turbulence develops further downstream. The Reynolds number at which the 
trails iti on occurs may be e^^ected to have the same order of magnitude 
as the Reynolds number for transition in flow along a flat plate. When 
the conditions are disturbed at entry the distance required for the velocity 
to take its final form is less but it depends on the amount of disturbance. 
When the flow is fully turbulent the inlet length 6 has been found to be 

8 - 1.3^6 a (3.15) 

which for R >10 is 8 > 13.86 which is considerably less than that for 
6 

laminar flow. Thus the error which results frem the assumption of steady 
flow in a pipe is less in the turbulent region than in the laminar region. 
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B) Flovr PMt a Flat Plate 
1 . Laalnar Flow 

Since the flow for some distance from the entry is similar to houndary 

layer flow past a flat plate a reasonable assumption is to treat the prob- 

1 2 

lem as boundary layer flow as was done by Carrier and Miles * . The drag 
force on a flat plate of width b and length i is 

D * 0.664 (3.16) 


From the dynamic analysis 

D = c^U^ a c^B = TimflRB 
The dan^jing constant n becomes 




1 . 328 (-|-) 

an 


1/2 J 


ITT 


(3.17) 


(3.18) 


The question which now must be asked is : What it the distance or 
length required for the boundary layer to disappear? Defining the bound- 
ary layer tWekness e as the distance for which u * 0,99 U then 

00 

(Schiicting,®pg. 122) 


e 



(3.19) 


Setting 6 equal to the radius of the pipe yields 

R 

6 - (^f)®- (3..20) 


which may or may not be a substantial portion of the fluid slug. Thus 
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in the ].aminar region neither the steady flow in a pipe approach or tbo 
■boundary layer approach may be a good approximation since neif icr is 
valid for the entire length of the fluid. 

2 , Turbulent Flow 

The drag force on a flat plate of ^r'.dth b and length I when the 

8 

boxindary layer is turbulent is (see Schlicting, pg. 536) 

« - 1/5 

D « 0.037 p u;b £(-^) (3.21) 

V 

Ecpiating this to *he viscous drag force 

D = c.U = c,R 0 
d * d 

one obtains 

, . o.arM^)'^^|)3^-V5,„.l|W5 ( 3 . 2 ,) 

This equation is valid for 

5 X 10^ < (|j)Rg < 10*^ (3.23) 

Spin Synchronous Mode 

In this mode the velocity of the fluid is not ccnsteuxt but oscil 
tory with respect to the ring. An approeush suggested by Leibold^ to ob- 
tain a dazoping constant is to \ise the results of Bbnita and Kbval^^, who 
analyzed the nutation dashing of a satellite with a coopletely filled 
viscous ring damper mounted on a plane parallel to the spin axis . They 
modelled the motion of the fluid as a fluid in an infinite pipe with the 
pipe executing nazmonlc motion, and then obtained the energy dissipation 

rate which leads to the dampl.ig constant. 
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To apply their aralysis to this problem it is assumed that the 
analysis is valid for a finite length of fluid, the energy dissipation 
rate is then averaged over m cycles to determine the average rate of 
energy dissipation. This average rate is then used to calculate the 
danping constant. From Bhuta and Koval the energy dissipation/unit- 
length at the end of the m^^ cycle is 


/sr, 

= 2U^“ ’^n 


s - , 4 y 2s2 
n=l (r^+C ) 


2 

-Ur mw 

exp I — — I -1 






(3.2U) 


where and I. are modified Bessel functions of the first kind and 
0 1 


I. = 

J 0 


I. = I . (>TiT ) 

J tJ 


(3.25) 


t - sa 

^ ~ V 


U is the maximum velocity of the tube with respect to the fluid, a is the 
radius of the tube, s is the excitation frequency and the r^ are the 
zeros of J^(r^), Letting g(0 "be the average amount of energy dissipa- 

p 

tion/cycle/unitmass/vel." one obtains for the damping constant 




(U 


(3.26) 


where 


I i.iA „ “ / 

;K ) = ?i| :::• 9 1 1 


2,2 


exp - 


Ur^mn^ 

“T 


-1 


,(3.27) 



g(^ ) is plotted for several values of m in Figure l6 . Since g(^ ) 
for m > 20 and g(^ ) for m = « are approximately equal, it is reasonable 
to just use the value of g(^ ) when m = » which is 


) = 


.1 y} 

' ^ 0^0 ^ 


( 3 . 28 ) 


The excitation frequency s 


IS 



0*01 



U. Test Data Analysis 


In October 1972 a series of tests were run at NASA/GSFC on the 
Helios damper. A total of 36 tests were run with four inertia ratios 
(0.337; 0.50, 0.51, 1.126) and two daji^jer locations (30 in. and 9 
above the satellite center of mass). In all of the tests the dan^jer was 
offset 0.25 in. from the spin axis. The test results and parameters 
as reported by Hraster^ are reproduced in Tables U.l and h. 2 . 

In each series of tests the satellite was balanced with the empty 
ring attached. Therefore, daring the tests one coxild consider that 
gravity is acting only on the slTig. With this assumption the effect 
of gravity on the nutational behavior of the satellite was determined 
in Section 2.3. 

In analyzing the test data the first thing which must be determined 
is in which mode, the natation synchronous or spin synchronous, the sat- 
ellite is operating. Because of the offset of the damper axis it is 
reasonable to assume the fluid is behaving like a rigid slug. In all 
of the tests the nutation angle time history appears to have an expon- 
ential behavior, hence it is assumed the satellite was in the spin 
synchronous mode. Using the development of Section 2 the effect of 
gravity on the time constant was removed and using Equation (2.30) a 
value of the damping constant t] and the transition angle 9 ^ were 
calculated. These resu' bs are given in Table U.3. Note that the time 
constant for zero gravity decreases with spin speed but the time constant 
for the tests increased with spin speed. Those tests for which there is 
no entry in the column there was no value of tj which would give 
time constant. However, for all of those except the last series of tests 
a 10^ change in the time constant would give- a reasonable value of t). 
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■><1 


m = 0.152 kg. 

R = 29 cm. a = 0.28 cm. 

= 0.25 7 = ^2 


V = 1.17 X 10 “^cm^/sec 

SEft. NO. 

A(slug-ft^) 

a 

l 3 =h/R 

i 

e 

0 . 50 - 4.1 

59.37 

0.337 

2.63 

6.35 X 10^ 

5.00-8.1 

65.85 

0.500 

2.63 

5.72 X 10“^ 

9.00-12.0 

65.85 

0.510 

0.788 

1 

-4 

5.72 X 10 

17.00-20.1 

40.05 

1.126 

1 

0.788 

. -4 

9.42 X 10 












ssq.KO 



c/o 

31.1 

76.7 

5.0 

61.6. 

383 

5.1 

60.9 

388 

5.2a 

60.4 

376 

5.2b 

60.4 

296 

5.2c 

60.4 

295 

6.0 

85.5 

554 

6.1 

81.2 

574 

7.0 

101.8 

816 

7.1 

102.2 

624 

8.0 

119.5 

1256 

8.1 

119.5 

803 

c/o 

31.5 

118 

9.0 

62.7 

379 

9.1 

61.7 

834 

10.0 

82.1 

1302 

11.0 

102 .1 

1468 

11.1 

99.7 

1390 

12.0 

120.0 

2032 

17.0 

62.3 

155.9 

17.1 

62.3 

191.6 

18 .o' 

81.1 

19.1 


i8.i 

19.0 

19.1 

20.0 

20.1 


81.2 

103.1 

102.6 

120.0 

120.8 


t^(g=o) 


14 

25 

14 


12.700 
15,200 

7,520 

9,300 

4,920 

7.700 


1,995 

1,970 

1,550 

1,550 

l,»+90 

1,680 

1,720 

1,310 

2,220 

1,420 


30,900 

20,500 

20.500 

13,700 
.8,780 
8,800 

8.500 





<? 


106.7 

- 

136.9 

4.1 

12.6 

7.8 

l4.l 

9.2 

6.54 

7.6 

6.89 

16.6 

3.96 

12:9 

3.90 

26.6 

2.83 

15.2 

2.78 


34.3 

5.4 

5.01 

5.2 

5.14 

5.3 

5.24 

- 

5.24 

- 

5.24 

6.6 

2.7 

5.7 

2.93 

4.9 

2.11 

6.9 

. 2.10 

3.4 

1.77 

5.6 

1.77 

3.9 

262. 

6.9 

23.4 

6.9 

24.6 

12.9 

10.5 

- 

5.99 

- 

6.34 

24.1 

4.19 

* 

3.24 


3.24 


2.16 


2.16 

am 

1.66 

• 

1.66 

- 

1.46 

m 

1.46 
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Using the analysis of Section 3 a value of has heen determined 
for each test and the calculations are presented in Table 4.4. Table 
4,4 shows that the value of T| calculated from the data is 4 to 5 times 
larger than the predicted value. It was originally thought that this 
difference was probably due to the offset but investigation has shown 
that the offset has very little effect on the time constant. Several 
possible reasons for this difference between the predicted and actual 
values of t) are; l) the method of calculating in Section 3 is not valid, 
2) the effect of gravity on the tests, 3) the mathematical model of the 
fluid behaving as a rigid sl.ug in the spin synchronous mocie, or 4) some 
combination of the above. 

In the fburth series of tests no value of the damping constant could 
be evaluated. It was originally thought that this was due to the offset 
since the offset angle is 1.5 deg. end the nutation angle was less than 
0.1 degrees. For nutation angles this small the offset did cause a 
decrease of about 25^ in the time constant but this is not enough to 
e^qplain the test results. Another contributing factor is that the 
nutation angles were so small that measurement of the time constant 
was difficult. 

For the first two series of tests the time constant has been scaled 
to the Helios satellite and is given in Table 4.5. A value of t) » 0.174 
was obtained by averaging the test data for 0 = 95 rpm. The correspond- 
ing time constant was then calculated. 

Since there is no test data for motion in the nutation synchronous 
mode there can be no comparison between actual and predicted time con- 


stants for that mode 
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I^(SLUG-FT^) 

Ig/lj = 

h(lN) 

0)(K£M) 

H (da25)ing const) 

t (flight condition) sec 
c 


0.385 

30.23 

95 

0.174 


59.37 

0.337 

30.0 

95 

0.174 


65.85 

0.50 

30.0 

95 

0.17 


48 , 500 . 


10,350. 


1 , 542 .. 












5. Summary 

A ring which is partially filled with a viscous fluid has been 

analyzed as a nutation damper for a spinning satellite. Since it was 

shown by Carrier and Mile^^^at the fluid behaves as a rigid slug for 

very small nutation angles the fluid haa been modelled as a rigid slug 

of finite length resisted by a linear viscous force . With these assunp- 

tions it has been shown that there axe two distinct modes of motion, the 

3 ^ 

nutation synchronous mode and the spin synchronous mode. For the 
symmetric satellite in the spin ssmchronous mode the nutation angle 
exhibits exponential behavior plus a small oscillation with the expon- 
ential Portion given by 

6 = 0^e (5.1) 

where 

In the nutation synchronous mode the cosine of the nutation angle 
exhibits exponential behavior 

t/t 

CO80 = cosQ^e (5»3) 


cn 


erT)(cT-l) 


For small angles (5*3) beccmes 


O T 


cn 


(5.4) 


(5.5) 
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The transition angle hetvreen the two modes is given hy 

tane = (5.6) 

^ ho^ 

Conroarisons of t and t with "exact” time constants obtained from • 
cs cn 

numerical integration of the equations of motion are given in Figures 
3-12. The agreement is good. 

The dan^jer was then analyzed for the asymmetric case and it was 
found that the two modes still exist. For the spin synchronous mode 




C7T) 


cs 2o^t(<Tj^-l)(<T2“l)+’l J [_ 






2 

cos 6 + 
o 




(5.7) 


In the spin synchronotis mode the slug oscillates in the tube while 

iDoving slowly around the tube. In (5.7) is the position abou.t which 

the slug is oscillating. Since P changes slowly with time t is a 

o cs 

slowly varying function of time, it oscillates between the two values 

of T obtained by setting P = 0 and P = jt/2 . For a design criteria 
CS o o 

one should use the maximum value of t . 

cs 

In the nutation synchronous mode for the asymmetric satellite we 

have 

0 m 
T 

cn 


lo^+crg-sl+lcg-c^lcosaXT 1/2 


(5.8) 


2o^8ign(ff^+ag-2) 

€?Ti(|c^+a2-2|+ja2-(^|) 


(5.9) 


The eq;uation for 6 is valid only for small nutation angles whereas 
the results for the nutation synchronous mode in the sysEoetrlc case 
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are valid for all nutation angles. Comparisons of and given 
"by (5*7) and (5.9) vith "exact" time constants are given in Figures l4 
and 15 . 

The effect of an offset of the center of the ring from the spin 

axis was investigated and was found to have only a very small effect 

on T and t , 
cs cn 

For a symmetric satellite an investigation was made of the effect 
of a stop in the tube. Since the behavior of the fltiid wL?n it en- 
counters a step in the tube is not known a very simple mathematical 
model wus used. The results show that the stop increases the amount 
of energy dissipation but no analytical result was obtaine 1 ;o predict 
the increase in energy dissipation. Some results are given in Sec- 
tion 2.4, 

Since t and t are a function of the damping constant t \ a meth'^d 
cn cs 

of calculating n is needed. In Section 3 several methods of calculat- 
ing T| are developed from a consideration cf the fluid dynamics. 

Analysis of the test results obtained from tests performed at 
NASA/jSFC on the Hellos satellite is given in Section 4. Before analyz- 
ing the results it was necessary to determine the effect of gravity on 
the behavior of the system. For the symmetric satellite it was found 
that gravity does not effect t but 

T {too gravity) 

T (test) - (5.10) 

cs 


where 


G - 




(5.11) 
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g is the ratio of the gravitational force to the centrifugal force 
which is the inverse of the Froude number, 

g - g/RO^ (5.12) 

Thus gravity can have a substantial effect on the test results. 

The test results show that during the tests all motion was in uhe 
spin synchronous node. Thus no comparison can be made with the theore- 
tical results for the nutation synchronous mode. The theoretical damp- 
ing constant developed in Section 3 was off by a factor of approxi- 
mately 4 or 5 from the value of ti calcnilated from the test resul.ts. 
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Appendix A 
Equations of Motion 


The system is assumed to consist of an asymmetric rigid body (the 
sateU-ite) and a circular tube of radius R which is attached to the 
rigid body at a distance h along the spin axis from the center 

of mass. The center of the ring is offset a distance 6 from the 
spin axis . Moving in the tube is a rigid slug which fi^JLs a portion 
of the tube, the fraction fill being 7, The only other assumptions 
are: l) the center of mass of the system and the center of mass of the 
satellite are coincident, 2) the motion of the slug is resisted by a 
linear viscous force, and 3) gravity acts only on the slug (the sat- 
ellite is statically balanced). 

Referring to Figure 1 the x, y and z axes are principal axes of 
the satellite and z is the spin axis. The u, v, z sys^em rotates 
about the z-axis relative to the x, y, z system such that the u 
axis passes throiogh the center of mass of the slug. Using the u, v, 
z coordin'^’oe system the equations of motion are obtained by equatirp 
the time rate of change of the angular momentum to the external moments 
and using Lagrange’s equation for the motion of the rigid slug in the 
tube. The angular momentum of the system about the satellite center 
of mass is 


H 


= cos 2 P + 1 * lo) - [-^^^3in20 + I* ]0) -I* (oH^Ae 

1 2 2 uu' u 2 uv* V uz z j Ml 

^ i 2p + 1* ] +[iA^_iA^ ^ 

12 uv^ u 2 2 w V vz z j ~v _ ' 


+ {-I oo-I oi+Coi+I (o-+^))e 
uz u vz V z zz z ~z 
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where A, B and C are the principal moments of inertia of the sat- 
ellite and the moments and nroducts of inertia of 

Cip 

the slug about the satellite center of mass . The I „ are given 

ap 

■in i^pendix B. 

The gravitational force is 


F = -mg N = -mg[sin0sin(^r+P)e +sin0cos )e +cos0e ] (A2) 


where 0, ^ and (ji are the Euler- angles of the satellite. The radius 
vector to the slug center of mass is 




[Rk+6cos (P-v)]^ 


6sin(p-v)e + he 


(A3) 


where 


k = 




(a4) 


The moment due to gravity becomes 
M = -mg 


^[-5sin(P-v 


)cos0-h sin0coc(\lr+p)]e 

•Ml 


+ [h sin0sin(ilr+P)-(Rk+5cos(P-v))cos6]e , 


(A5) 


+[ (Rk+6cos (P-v ) )sin0cos (t+P )+6sin(p-v )sin0sin(^lr4^ )]e 

<vZ 




n 

I ^ 


The kinetic energy of the fluid slug is 


■K 


TL, „ = —mv^ + ‘CO = + ^mv'^ +mv -(coxr ) (a6) 

2 c 2 'v '\FS 2 cr ~cr ''V' '^<c 


1_„2 


where v is the velocity of the center of mass of the slug, v^^ 

" <vcr 

is the velocity of the center of mass of the slug relative to a coordinate 


i 




i 
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system whose origin is at the satellite center of mass and whose axes 
are parallel to the u, v, z axes, (u is the angular velocity of the 
fltiid slug and is the angular momentum of the slug. Using (A3) 


V = -8^[sin(6-v)e +cos(p-v)e J 


(AT) 


The kinetic energy becomes 




I - * 2 * 2 */ •^2 * * / *. * / 

CO +I CO +I (co +0) -21 CO CO -21 co (co 4^)-2I co (co "t^)] 
2 UU U W V ZZ Z UV U V UZ U Z VZ V z 


-m6p[h(co sin(P-v)-co cos(P-v))+Rk(co +0 )cos (6-v)+6(co +6/2)] 
V u z z 


(A8) 


The potential energy is 

V = mg r^*N 

^ ~c <vz 

V = mg Kk sin9sin(Tlf+6)+Bsin0sin(ilf-v)+h cos0] (A9) 

The generalized force due to the line i.scaus force is 

Op = -c^rS* (AlO) 


In the development of the approximate solutions it is advantageous to 
use dimensionless variables and constants. The angular rates and time 
are made dimensionless using the initial spin rate. Let 


T = Qt 


CD ® 

V 

£ 0 ^ = flr 


(All) 


(A12) 
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A suitable set of (dimensionless parameters are 


= C/a 

cTg = C/B 

B = h/R 
8 = 8/R 

Tj = 

mR^ 

e =: —— 

A7 


g 


= -S_ 
Rn^ 


(A13) 


In addition to these we have 7, the fraction fill, and v. 6 is a 
small parameter which is the ratio of the moment of inertia of the 
tube filled with fluid to one of the transverse moments of inertia 
of ‘the satellite. It was chosen in this manner so that e would 
remain constant when varying 7 or 7 . ti is a dimensionless dairping 
parameter, g measures the effect of gravity and is actually the inverse 
of the Froude number. Substituting (A11)-(A13) into (AI), (A5), (a 8), 
(A9) and (AlO) gives 
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+ fl P-I <1-1 
‘ mr uv^ uv' 


')i]s. 


(1« ) (1-P15) 


r (l-a ) 

+ |. [ ^sin2p]p + [ 

+ f<r,r + [-I _P-I a+I (r^e')n e 


COS 23lq 


uz vz zz 


= V»2 


O' - ^ 


= -€7g ([-6sin(P-v)cos0-t) sin0cos(t^)l^ 


[■b sin0s±n(t-t^ )- (k+&cos (P-v ) )cos9]e^ 


+[ (k+6cos 0-v ) )sin0cos (t-<# )+5sin(P-v )sin0sin(t-*^ )]^l 

^ = kl P^+i q^+i 0^-21 ^q-2i ^(r4P')-2i <l(r^’)] 
p 9 p*- iitr w zz u\r ur w* 


„2«2 2'- Mvr zz 

fl 


-8P'['b(q sin(P-v)-P cos(P-v))+k(r+P* )cos(p-v)+6(r+P'/2)] 


-X_ = g[k sinesin(iHP)+ sin9sin(t-v)+'b cosS] 
jiirV 


A_ = 

n D r 


T,2ft2 

nR fl 



Applying 
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dH 

dt ^ ~ 


and 


^ St 

dt ^ 


S-| 


and noting that 


St St 
3F = 5p 


p,q = const 


* -^55 


since p and q depend inq)licitly on p via 


0 ) = CO cosP + CO sin3 

u X y 


(0 = -CO sin3 + CO cosP 

V X y 


the equations of motion become 


(1+CT,^) 






7 ^cos 2 P+I ]p» + [- sln2g-I ]q’ 

2 uu ^ ''2 


+ P^(r+e«)p' + (r^*)S^ - qflj, + \ 


(A21) 


(A22) 


(A23) 


(A24) 


(A25) 
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[- — ^sin2g-I^^]p' + ^cos20+I^k’ 


I_(r'4p") = [-5i^cos2p+rjpP* + [-I^-(l-a^)sin2p]q?’ (A26) 


vz 


uv 


+ r^(iH^')P’ + pfljj - + 1\ 


•I p’ - I q‘ + (o,+i )r' + I 3” * I’ pe' + I’ qP* 
uzr vz^ zz zz uzr vz 


(A27) 


+ qH^ - PS^ + 


-bl^' + (l+Bk cos(3-v))r* + p" = - g k sinPcos (\lf4^) 


+ (I -1 -r )pq - bk <i(r4^’) (A28) 

^ uu w uv 




- “a 


«a = 


Afi^ 


I» 

ccu) 


bl 

Ofco 

■IT 


Some of the terms in (A28) have been simplified by substituting 
for the their values given in ^pendix B« 

QU) 

The relationship between the ihiler angle rates and p, q, r is 
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p = 0cos(^) + <jisin0sin(tt^) 

q = -esin(tt^) + <l>sin0cos (\|r+p ) 


= ijr + ^cos® 


(A29) 


IH 


Appendix B 
Moments of Inertia 


vcu 


w 


zz 


uv 


uz 


vz 


uu 


w 


zz 


uv 




vz 


oe 


+ B^sin^(P-v)] 

2 

m[ ( toASZ l + 26R k cos(B-v)+B^cos^(p-v)] 

m[R‘"+25R k cos(P-v)+S^] 

-mBsin (p -V )[R k +6cos (p -v )] 
nih[R k+Bcos (P-v)] 

-mh Bsin(P-v) 

(i.s ^)/2 + ^2 6^sin(P-v) 

( H . g i ^7 )/2 + + 2§k cos(P-v) + 5®cos^(p-v) 

1 + 26k cos (P-v) + 6^ 

-6sin(P-v)[k+8cos(p-v)] 

l3(k+6cos(p-v)] 

-li6sin(P-v) 
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■XJf 




1 

zz 



€^sin2(0-v) 

-2Sk sin(p-v) - 8^sin2(p-v) 
-2Sk sin(p-v) 

-k5cos(p-v) - S^cos2(p-v) 
-IsSsinCP-v) 

-T35cos(P-v) 
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